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ABSTRACT 



We present a new diagnostic allowing one to measure the anisotropy of ion 
temperatures and non-thermal velocities as well as Doppler shifts with respect 
to the ambient magnetic field. This method provides new results, as well as 
independent test for previous measurements obtained with other techniques. Our 
spectral data come from observations of a low latitude, on-disk coronal hole. A 
potential field source surface model was used to calculate the angle between the 
magnetic field lines and the line of sight for each spatial bin of the observation. 
A fit was performed to determine the line widths and Doppler shifts parallel and 
perpendicular to the magnetic field. For each line width component we derived 
ion temperatures T^± and and non-thermal velocities v nt> ± and u n t,||- ^1,± 
was cooler than off-limb polar coronal hole measurements. T^y is consistent with 
a uniform temperature of 1.8 ± 0.2 x 10 6 K for each ion. Since parallel ion 
heating is expected to be weak, this ion temperature should reflect the proton 
temperature. A comparison between our results and others implies a large proton 
temperature gradient around 1.02 R Q . The non-thermal velocities are thought 
to be proportional to the amplitudes of various waves. Our results for v nt; ± agree 
with Alfven wave amplitudes inferred from off-limb polar coronal hole line width 
measurements. Our v n t t \\ results are consistent with slow magnetosonic wave 
amplitudes inferred from Fourier analysis of time varying intensity fluctuations. 
Doppler shift measurements yield outflows of ~ 5 kms -1 for ions formed over a 
broad temperature range. This differs from other studies which found a strong 
Doppler shift dependence on formation temperature. 
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Introduction 



Wave- and turbulence- driven models of coronal heating and solar wind acceleration 
propose that the necessary energy is car ried to larger h eights via plasma waves driven by 



agitation in and below the photosphere (jCranmer 



20021 ). Attempts to theoretically describe 



the dissipation of these waves and how that energy is converted into particle heating 
represent a major area of research in solar physics. Observations of the temperature, wave 
properties, and flow velocities of the plasma can constrain these models. In the corona, 
where solar wind acceleration begins, emission line profiles can be used to infer many of the 
initial properties of the solar wind. Further out, in-situ measurements near 1 AU provide a 
detailed description of the plasma. 

The bulk motions induced by waves are expected to be anisotropic with respect to the 
ambient magnetic field. The two main types of wave in the solar corona are Alfven waves 
and slow magnetosonic waves. Alfven waves are transverse waves and can be observed 
spectroscopically through non-thermal line broadening perpendicular to the magnetic field. 
here is strong evidence that Alfven wave s exist throughout the Sun from the chromosphere 



( De P ontieu et a 



( iBelcher fc Davis 



2007) into the corona ( ITomczyk et al. 



20071 ) and into the solar wind 



19711 ) . Spect roscopic measurements indicate that 



these waves can be 



important for coronal heating ( iHahn et al. 



2012 



Bemporad fc Abbo 



20121 ). 



On the other hand, slow magnetosonic waves cause non-thermal broadening in the 
direction parallel to the magnetic field. These are essentially sound waves, but they are 
modified in the solar corona by the strong magnetic pressure which constrains the fluid 
motion of the wave to be nearly along the magnetic field. Because they are compressive, 
these waves may also be observed as time varying line intensity oscillations. Intensity 
fluctuations consis tent with slow magn etosconic waves have been reported in a number of 



observations (e.g.. iBanerjee et al.ll201ll ). 
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It is also known that the ion temperature can be anisotropic in the solar wind. For 
example, in-situ measurements near 1 AU indicate that proton temperatures are anisotropic 



with T p j_ > T p ii (IMarsch et al. 



2004 ). Such anisotropic hea ting can be caused by cyclotron 



resonance with high frequency waves (ICranmer et al. 



1999T ) o r by stochastic 



large amplitude fluctuations transverse to the magnetic field ( jChandran 



reating from 



2010h 



However, closer to the Sun the ability to detect ion temperature anisotropy and 
distinguish between Alfven and magnetosonic waves through emission line profiles has been 
limited because spectroscopic observations usually only detect one component of the line 
width. Typical line width measurements in off-limb spectra of coronal holes or quiet Sun 
regions are primarily sensitive to perpendicular broadening. This is because the line of sight 
is nearly perpendicular to the magnetic field lines at the point where the line of sight passes 
closest to the Sun. Since the electro n density n e . drops exponentially with height, and line 



intensity is proportional to (e.g. 



Hahn et al 



20101 ). any varations of the angle along 



the line of sight has little effect on the measured width. Thus, such measurements provide 
a useful diagnostic for Alfven waves and perpendicular ion temperatures, but provide no 
information about sonic waves or the parallel ion temperature. 



I 

Several attempts have been made to overcome these limitations. iKohl et al. 



(1998) 



observe d a polar coronal hole above 1.5 R & with the Ultraviolet Coronagraph Spectrometer 



(UVCS; 



Kohl et al 



19951 ). By combining line width and Dopper dimming diagnostics for an 



O VI line they were able to measure the perpendicular ion temperature T^± as well as put 
upper and lower bounds on Ty. They determined that 7]^ > Ty down to at least 2.2 R & . 

An alternative method to separate parallel and perpendicular broadening has been to 
measure the variation of line widths from the center of the solar disk to the limb. The 
principle of this method is that, on average, the line of sight looks down the field lines of 
large scale magnetic loops at the center of the disk, while near the limb the observations tend 
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to look across the field lines. The results from such studies have been mixed. Feldman et al. 



by 



(119761) found no change i n the line width, whereas center -to-limb broadening has been found 



Roussel-Dupre et al. 



(119791 ) and 



Erdelyi et al. 



(Il998f ). Even if a clear trend could be 



discerned, the precise separation of the parallel and perpendicular broadening would remain 
ambiguous because the magnetic field direction is not known except in an average sense. 

Flow velocities provide another constraint on solar wind models that complement the 
temperature and non-thermal velocity data. The flow of material from the Sun into the 
solar wind can be inferred from Doppler shifts. The ability to make precise measurements of 
the Doppler shift depends on the precision with which the rest wavelength on the detector 
is known for the observed spectral line. Since most spectrometers do not carry an onboard 
calibration lamp it is usually not possible to use the laboratory wavelength for this. One 
way to determine the rest w avelength is to study the center-to-limb variation of the line 



center (IPeter fc Judge 



19991 ). This method assumes that at the limb any non-perpendicular 
motions cancel out when summed along the line of sight. The disadvantage of this method 
is that the Doppler shift is measured only in an average sense. Another method is to assume 
that there are no flows in the chromosphere and measure all wavelengths with respect to a 
chromospheric line. Since there probably are flows in the chromosphere this method is only 
accurate when the velocity in the corona is large relative to any chromospheric velocity. 

Here, we describe a new diagnostic which allows one to measure both the parallel and 
perpendicular components of a line width and also the flo w velocity. We use an observation 



Culhane et al. 



20071 ) taken of a 



with the Extreme ultraviolet Imaging Spectrometer (EIS; 
coronal hole at a relatively low latitude. The slit spanned a region where there was a large 
variation in the angle between the line of sight and the magnetic field. A Potential Field 
Source Surface (PFSS) model was used to trace the magnetic field lines passing through 
each spatial bin of the EIS data. As we explain later, this allowed us to calculate both the 



-6 - 



line width and Doppler shift as a function of the inclination angle between the line of sight 
and the magnetic field and thereby separate their parallel and perpendicular components. 
For the line widths, we analyzed the two components separately to estimate the partition 
between thermal and non-thermal broadening. From the Doppler shifts we determined 
the bulk flow by assuming that there is no motion perpendicular to the field lines. As we 
explain below, this allows us to fit for the rest wavelength and thereby infer the parallel 
velocity from the Doppler shift relative to that wavelength. 

The rest of this paper is organized as follows: Section [2] describes our analysis including 
details of the observation, fitting of spectral lines to determine the line widths, the PFSS 
magnetic field model, and separation of the parallel and perpendicular components. Results 
are presented in Section [3j In Section [3TT1 we characterize the observed region in terms of the 
electron temperature and density. Sections 13.21 and 13.31 then discuss the ion temperatures 
and non-thermal velocities extracted from the perpendicular and parallel line widths and 
compare these results to measurements using other methods. Our Doppler shift results are 
presented in Section 13.41 Section H] summarizes our conclusions. 



2. Analysis 



2.1. Observation 



We analyzed archival EIS data of a low latitude coronal hole. The observation was 
performed on 2007 September 26 at 14:17 UT. The 1" slit was rastered across 36 positions, 
each having an exposure time of 90 s. The data come from a 456" long portion of the 
EIS slit. The observation was centered vertically at -537" and the slit was rastered from 
X = 321" to 349". Figure [T] shows the location of this o bervation superimposed on an image 
from the Extreme ultraviolet Imaging Telescope (EIT; iDelaboudiniere et al.lll995l ). The 
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image was taken at about the same time as the EIS data. The coronal hole we observed 
appears to be part of the South polar coronal hole. 

The data were prepared using standard EIS procedures to remove spikes, warm pixels, 
and CCD dark current, and calibrate the intensity scale. S ystematic dri f ts in t he wavelength 



Kamio et al. 



3). 



(120101 ). In order to 



scale were then corrected using the method developed by 
increase statistical accuracy, we binned the 1" square pixels into bins 9" horizontally and 
19" vertically. The size was chosen so that each bin would contain the same number of 
pixels. Thus, the effective exposure time for each bin was about four hours. Howey er this 



is still rather low compar ed to the detailed line width studies by 



Hahn et al 



tom and 



Bemporad fe Abbol (120121 ). which each used data binned to effective exposure times of over 
twenty hours. 

Figure |2] illustrates the binning as a grid over an unbinned EIS image in the Fe XII 195 A 
line. The observation covers several structures, including the coronal hole, the quiet Sun, 
and several bright points. As we are interested only in the coronal hole bins we have 
excluded from our analysis bins that appear to include bright points and quiet Sun material. 
The remaining bins selected for analysis are marked by the dots in Figure |2j 



2.2. Line Widths 



Gaussian profiles were fit to each spectral line in order to determine the line width AA. 
In some observations, particularly of active regions, it has been s hown that a Dop pler-shifted 



component c an cause an apparent broadenin g of the line profile (ITian et al. 



201lh . For quiet 



Sun regions, iMcIntosh fe De Pontieul (120091 ) found blue-shifted components of an Ne VIII 
having a velocity of ~ 50 kms -1 with a relative intensity of ~ 5%. The shifted components 
were strongest in the supergranular network. We tested our line profiles for asymmetries by 
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comparing the integrated intensity on each side of the line center determined by a single 
Gaussian fit and did not find any significant asymmetries. One possible reason for this is 
that the broad spatial binning covers both network and internetwork regions, which washes 
out the asymmetry. 

Some studies have also shown that line profiles observed on disk have broadened w ings 



and requi re double G aussian fits (IKjeldseth Moe & Nicolas 



However, 



1977 



Wilhelm et al. 



20071 ) 



Peterl (120011 ) measured the relative intensity of the broad and narrow components 



as a function of line formation temperature and showed that the intensity of the broad 
component goes to zero at coronal temperatures. We therefore expect single Gaussian fits 
to be sufficient for our study which focusses on coronal lines. We have checked this by 
comparing double to single Gaussian fits for some lines, and found no improvement in the 



fits. Other studies of EIS data from coronal holes have a 
well approximated by a single Gaussian (e.g., iTian et al- 



so fo und that the line profiles are 



20101 ). 



Table Q] lists all of the lines used in our analysis; the transitions selected for the line 
width study are marked with an asterisk. We have avoided using blended lines except for 
the Mg VI 270.40 A, Fe x 256.27 A, and Fe xn 195.12 A doublets, which are self-blends. 
In order to accurately measure AA, we fit these lines with a double Gaussian constrained 
so that both components have the same width and a fixed separation between the two 
centroids given by the known wavelengths. 

The spectral line full width at half maximum AAfwhm depends on ins trumental 



broadening AA; nst , the ion temperature T[. and the non-thermal velocity v ut (Phillips et al 
20081 ) as 



AA 



FWHM 



AAL t + 41n(2) 



A 



2k B Ti 2 



1/2 



where M is the ion mass, k-Q is the Boltzmann constant, A is the observed wavelength, and 
c is the speed of light. The instrumental width is known to vary along the length of the slit. 
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Yound (120 



Hara et al 



a ) m easured this variation by studying quiet regions near the solar equator. 



( 2011 ) obtained similar results by comparing the widths of the Fe xiv 264.78 A 



the Fe XIV 5303 A line. Here, we 



Young (j2011al ). The remaining line 



line observed by EIS with ground based observations o: 
subtract the instrumental widths using the values from 
width is the sum of the thermal and non-thermal velocities. This width can be expressed in 
terms of an effective velocity, defined as 



VeS 



M 



(2) 



2.3. Magnetic Field 



A PFSS model was used to trace the magnetic field lines in 



;he coronal hole 



( Schatten et al. 



1969 



Wang fe Sheeley Jr. 



1992 



Schrijver fe Derosa 



20031 ). The model 



assumes that there are no currents between r = R & and the source surface at r = R s . The 
radius of the source surface is arbitrary, but a value of -R 9 = 2 . 5 has usually been found 



to give the best agreement with observations (IHoeksema et al 



1983). In the region between 



Rq and R s the magnetic field B can be determined by solving the Laplace equation. The 
lower boundary condition is to require that the solution at r = R Q matches measured solar 
photospheric magnetogram data. The boundary condition at r = R s is that the field lines 
become radial. 

Although the PFSS model makes significant simplifications, it generally does a good 
job of reprod ucing the large scale structure of the solar magnet i c field, particularly for 



coronal ho 



Riley et al 



es (TWang fe Sheeley Jr. 



1992 



Neugebauer et al. 



1998 



Schrijver fe Derosa 



2003 



20061 ). This has been tested by tracing the open field lines, defined as those 



that reach R s , back to the solar surface and comparing these to the observed coronal hole 
boundaries. The coronal hole boundaries predicted by PFSS and full MHD models are in 
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good agreement with observations. This suggests that PFSS models provide a reasonable 
approximation for coronal hole magnetic fields. 

Here, we used the model implemented by th e PFSS package available in solarsof\\ and 



described in detail by 



Schrijver fc Derosa! (120031 ). The data for this model is available at 



six hour intervals. We chose the one at 2007-09-26 12:04, which was closest in time to our 
observation. An inspection of PFSS results from before and after this time showed that the 
magnetic field lines in the coronal hole appeared to be stable during the observation. In the 
analysis we corrected for the very small shift in viewing angle due to the rotation of the 
Sun between time of the PFSS model and the observation time. 

In order to find the angle between the magnetic field line and the line of sight for each 
spatial bin in our data we used the PFSS package to trace the magnetic field line passing 
through the center of each bin at a lower height R\. Then we found the coordinates of 
each of these field lines at a larger height i?2- The unit vector pointing between these two 
coordinates gives the average magnetic field direction b between R\ and i?2- We label the 
direction outward from the Sun parallel to the line of sight as the z direction. The angle 
between the magnetic field and the line of sight is then a = cos _1 (b • z). 

Figure [3] shows the field lines passing through the bin centers. The figure shows that a 
depends on the radii over which the angle is calculated. The field lines become straighter, 
and a becomes more constant, for larger heights. Here, we chose R± = 1.02 R Q and 
i?2 = 1-05 Rq for the analysis. These heights are based on our density and temperature 
measurements, which are consistent with observing at a height of about 1.02 R & with an 
emission measure scale height of ~ 0.03 R Q (see Section [3TTj) . We have also tested our 
results for line widths and Doppler shifts using height ranges of 1.02 - 1.10 R Q and 1.01 - 
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1.05 R & . The results for the 1.02 - 1.10 R & range compared to 1.02 - 1.05 R Q were nearly 
the same. For 1.01 - 1.05 R Q there were larger differences compared to 1.02 - 1.05 R & , but 
the changes were still generally within the la uncertainties. 

There are also small variations of the angle a in the plane of the observation with 
height. Due to the large spatial binning a varies from the nominal value determined at the 
bin center by ±3° on average. 

Figure [3] also demonstrates that a simpler assumption of a radial magnetic field would 
not be accurate. The field lines are radial only far from the boundaries of the coronal hole. 

2.4. Anisotropy 

The line width reflects the broadening due to v e g along the line of sight. The velocity 
of any particular ion in the observation can be described by v = vnh + v±q where b is the 
unit vector along the magnetic field and q _L b. Thus, the projection of v along the line of 
sight z is v los = v • z = wj| cos a + v± sin a. We expect v\\ and v± to be normally distributed 
with variances v 2 ff « and v^ s ± . Since v • z is a linear combination of v\\ and v± the variance 
along the line of sight v 2 ff comes from the convolution of the two Gaussians, which leads to: 

<4 = v cs,\\ cos 2 (a) + v* Sj± sin 2 (a). (3) 

This is the quantity that we measure from the line width. 

By fitting v e g versus a we can extract the components v e g t \\ and f e ff,_i_ from the measured 
v e ff. A least squares fit was performed separately for each ion species. An example of a fit 
for the Fe XI lines is shown in Figure HI As we discuss later, there are some indications that 
our coarse spatial binning covers different structures within the coronal hole within each 
bin. If this is the case then a longer observation might remove some of the scatter in the 
plot by permitting smaller bin sizes. Table [2] summarizes the results for all the ions. The 
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inferred t^gii and v e g t ± include both thermal and non-thermal broadening in the same way 
as described by equation [2J The evaluation of these contributions will be discussed in more 
detail in Sections 13.21 and 13.31 

3. Results and Discussion 

3.1. Density and temperature 

The electron density n e and temperature T e are useful for interpreting the line width 
results. The density averaged over the selected bins was determined from several line 
intensity ratios. From the Fe vill 185.21 A/186.60 A intensity ratio we inferred an 
average density of n e = 1.0 ± 0.2 x 10 8 cm" 3 ; the Fe IX 189.94 A/188.49 A intensity 
ratio gave n e = 2.1 ± 0.7 x 10 8 cm" 3 ; and the Fe xill 203.8 A/202.04 A ratio implied 
n e = 2.6 ±1.3 x 10 8 cm -3 . Here and throughout all uncertainties are given at an estimated 
la level. Note that the Fe xm 203.8 A line is a blend of four Fe xm lines, all of which were 
included in our analysis. The most important lines are the ones at 203.796 A and 203.827 A, 
but for n c < 5 x 10 8 cm -3 contributions from the lines at 203.772 A and 203.835 A are also 
important. 

These values are within a factor of a few, which is interesting since the Fe vm, IX 
and xill ions are formed over a range of temperatures. Both the transition region and 
the corona lie along the line of sight, so the observation looks into a temperature gradient 
with T e increasing with height. In contrast, the density should be decreasing with height. 
Here we find that the densities roughly agree, or possibly follow a temperature dependence 
opposite to what is expected, although the uncertainties are too large to say so definitively. 
This agreement suggests that the ions are formed in the same volume. However, another 
possibility is that there are unresolved hotter denser structures that skew the density 
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measurements. 



For an on-disk observation there will be material at a range of different temperatures 
along the line of sight. In order to measure T e we performed a differential emission measure 
(DEM) analysis. The DEM 0(T e ) describes the amount of material along the line of sight as 
a function of T e . In terms of 4>(T e ), the intensity of an emission line emitted by a transition 
from level j to level i is given by 



j' 



j- J G(T e )0(T e )dT e 



(4) 



where G(T e ) is the contribution function and describes the level populations, ionization 
balance, elemental abundanc e, and radiative decay rates. These data are available in the 



CHIANTI atomic database (iDere et al. 



1997 



Landi et al 



2012al ). Since we started this 



work the CHIANTI database has been updated from Version 7 to Version 7.1, but these 
updates have little effect on the plasma DEM and do not affect the analysis of line widths 



and Doppler shifts (ILandi 



20121 ). Given G(T C ) and a set of measured line intensities Iji, it 



is possible to invert equation (jlj) to find 4>{T C ). 



We calculated the DEM using the technique described in ILandi fc Landinil (119971 ). The 
lines used for this analysis are listed in Table [1] and include some ions formed at lower 
temperatures than those considered in the line width analysis. For the inversion it was 
necessa ry to set the behav ior of 0(T e ) at both the high and low temperature ends. As was 



done in 



Hahn et al 



(1201 ll ). at high temperatures we assumed that 0(T C ) = at logT e = 8 
(here and throughout all temperatures are reported in Kelvin). This is a reasonable 
condition since there should not be any material at higher temperatures in the observation. 
Due to the presence of the transition region and chromosphere, 0(T e ) will not be zero at 
low temperatures and it was not possible to anticipate the low temperature behavior of the 
DEM prior to performing the analysis. For the DEM reported here we have set <f)(T c ) at low 
temperatures to a constant value determined by the He II line, which is the coolest line in 
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our data. We also tested a variety of other conditions and found that for logT e > 5.6, the 
DEM does not vary significantly. This is because the contribution function for the coronal 
lines goes to zero at low temperatures, so that temperature range contributes nothing to 
the coronal line intensity integrals. As an additio nal check we also p e rform ed the DEM 



Hannah fc Kontarl ((2012). Both DEM 



analysis using the regularized inversion method of 
inversion methods agree to within their respective uncertainties for logT e > 5.6. Thus, 
despite the poor constraints at low temperatures the DEM appears to be valid for coronal 
temperatures. 

We calculated the DEM for various positions within the observation and found that 
the DEMs from all the selected coronal hole bins were similar. Figure [5] shows the DEM 
calculated for a typical position in the coronal hole, X = 325" and Y = —525". In the 



figure the solid line indicates the DEM calcul ated using the 



and the crosses represent the DEM from the 



Lancli fc Landinil (119971 ) method 



Hannah fc Kontarl (120121 ) method. The DEM 



shows that the observed emission comes from a broad range of temperatures, though there 
is a peak at logT e w 6. The dots in the figure correspond to the various lines used in 
the DEM analysis. The error bars on these points are from the uncertainty of intensity 
only, and neglect possible additional errors from atomic data and unknown blends. Their 
position on the temperature scale corresponds to the DEM-averaged temperature at which 
each line is emitted: 



logT t 



To determine the T+, we used the 



fG^(T e )0(T e ) log T e dT e 
f Gji(T )<f>lT )dTr 



(5) 



f 

Landi fc Landinil (119971 ) method DEM results. Since 



the temperature dependence of is mainly due to the ionization balance, the effective 
formation temperature T t is about the same for all lines from a given ion species. Table [2] 
lists T t for the ions used in the line width analysis. 



The DEM has a small peak at about logT e = 6.0 and drops off towards higher 
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temperatures, but at lower temperatures cj) does not drop due to presence of the transition 
region along the line of sight. The extension of the DEM up to log T P ^6.2 may be du e 
to the limited temperature resolving power of the DEM inversion. 



LandietaL 



(l2012bl ) 



studied the ability of a DEM technique to measure isothermal plasmas and found that 
the method could not resolv e two isothe r mal c omponents if they were separated by less 



than about AlogT e = 0.20. lHahn et al. 



(j201l[ ) use the same DEM technique as here to 



study an off-limb coronal hole observation, and found that the full width half maximum for 
the peak in 0(T C ) was about logT e = 0.15. This raises the possibility that the lines from 
higher char ge states, su c h as Fe XII and xm, include contributions from hotter structures. 



Previously, 



Hahn et al. 



( 1201 ll ) measured the DEM of a polar coronal hole and suggested 
that the high temperature tail could be due to quiet Sun coronal material intervening along 
the line of sight. However, this explanation does not apply to the present observation where 
the PFSS model shows that the closed field lines are bent away from the coronal hole and 
do not intersect the line of sight. In this case a more likely explanation is that the large 
spatial binning sums over temperature variation due to small scale structures within the 
coronal hole. 

Because the intensity is proportional to n\, the observed emission will be dominated by 
the points along the line of sight with the greatest density. The measured log T e ~ 6 implies 
a density scale height of ~ 0.06 R@. Since the scale height for nl is half the scale height for 
n e , we infer that the measured intensities come primarily from a height range of ~ 0.03 R Q , 
assuming a constant temperature. Our density measurements can be used to specify the 
height in the coro na w here the emission is produ ced. Off-limb density measurements by 



Hahn et al. 



(120101 ) and iBemporad fe Abbol (120121 ) showed that the measured density of 



n e R!2x 10 8 cm 3 corresponds to heights of about 1.02 R & to 1.05 R & . Thus, the observed 
emission appears to come from approximately these heights. 
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3.2. Perpendicular Broadening 



The component i> e ff,_i_ perpendicular to the magnetic field depends on T^± and v nty 
through a relation analogous to equation (j2J), namely 



VeS,J 



2fc B T iv 
M 



+ v 



nt,_L 



(6) 



In coronal holes T it j_ has been observed to be greater than T P a nd varies with th e ion 



charge to mass ratio (jEsser et al. 



1999 



Landi fc Cranmer 



2009 



Hahn et al. 



2010) • The 



perpendicular he ating may be due to such thin gs as ion cyclotron resonance with high 



frequency waves ( 



Cran 



( Voitenko k, Goossen 



mcr 



2004 : 



2002 



Hollwee 



Chandran 



201 



2008), or stochastic heating by turbulence 
The non-thermal velocity perpendicular 



to the magnetic field is thought to be 



(IDoyle et al. 



1998 



Banerjee et al. 



1998 



oroportional to the amplitude of Alfven waves 



20111). 



Upper and lower bounds on Tj i and an upp er bound for t> nt ,_L were found from t> eff 



using a method based on that of 



Tu et al 



(119981). Their approach assumes that all the 



emission comes from the same volume so that the fluid motion from the waves, seen as 
Vnt,±, is the same for all the ions. In our observation the volumes are large and could 
encompass different structures, which might make this assumption invalid. However, as 
discussed below, our results do appear to be consistent with t> nt j_ the same for all the ions. 



T he upper bound for T^± was determined by setting v nti ± = 0. To find the lower 



bound 



Tu et al. 



(119981 ) assumed that for the narrowest line the width was entirely due to 
non-thermal broadening and applied that value of f n t,_i_ to determine T- h ± for the remaining 
lines. We found that assumption leads to a lower bound on T^± > T e for most lines. Since 
we expe ct that the i o ns wi ll either be in equilibrium with the electrons or be heated, we 
followed 



Hahn et al 



(120101 ) and derived a tighter constraint on the lower bound of 2] j_ by 



assuming 7] j_ > T e . The maximum v nt ± is given by the minimum width after subtracting 
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this lower bound of T e from t> e fr,_i_- For this analysis we took logT e = 5.67, since this was the 
lowest value of logT t for these ions (Table [2]). 

Figure |6] shows the upper and lower bounds for T^± as a function of charge to mass 
ratio q/M. For low q/M ions, even the lower bounds for T- h ± are greater than T t , which 
we take as a measure of T e for each ion. The inferred T^± decreases for q/M > 0.16 . 



The se results can be c ompared to similar measurements by 



and 



Hahn et al 



Landi fe Cranmerl (120091 ) 



( 120101 ) from off-disk spectra. The shape of T^± versus q/M is similar, 
with a high temperature at low q/M that drops at higher q/M. However, in the present 
observation the ion temperature for the lowest q/M ions was smaller than that found in 
off-limb observations where log 7]^ was > 6.5. 

A possible explanation for these differences is that here, the data come from a lower 
height in the corona. The ions are expected to be cooled through collisions with the protons 
which are expected to be cooler than the ions. The density here is about 2 x 10 8 cm" " 
compa red to 8 x 10 7 cm -3 at the lowest height in the off-limb observation of 



Hahn et al. 



( 120101 ). Thus, collisional cooling is more important and this may make T{± smaller here. 



The coarse spatial binning of our data introduces a possible systematic error through 
the assumption that f n t,± is the same for all the ions because different structures might 
not have the same f n t,_i_- We can test the validity of this assumption by checking the 
consistency of T^j, for ions formed at different T t . For example, Si VI at \ogT t = 5.67 and 
Fe XI at logT t = 6.08 have nearly the same q/M and their T- h ± are in reasonable agreement. 
Another example is the similar level of agreement in found between Si VI and Fe xill. 
The consistency of the results suggests that our assumptions in the analysis were justified. 
However, we should note that the possible presence of different structures in the data is not 
ruled out. It is possible that there are multiple structures, but they happen to have the 
same v nt> ±, or that our measurements are insensitive to small inhomogeneities. 
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The above analysis for T{ t ± yielded an upper bound for u n t,i of 2 4.9 kms 1 
undamped Alfven waves on open field lines we expect v nt ± oc n^ 1 ^ (Moran 



2001|). 



or 



Off-limb observations have shown that v nt ± follows this dependence, which supports the 



interpretation of v n t ± as Alfven waves. Those obse r vations have a. 



are undamped below ~ 1 . 15 R & (iDoyle et al 



1998 



Bemporad fc Abbo 



Banerjee et al. 



so shown that the waves 



1998 



Hahn et al 



2012 



20121 ). Thus, we can use this proportionality to scale our measured 
Vvt,± for comparison with off-limb observations. For n e w 2 x 10 8 cm -3 here and a typical 
value of 8 x 10 7 cm -3 at 1.05 R Q in off-limb observations, the scaling factor is ~ 1.26 and 
the expected upper bound for i at 1.05 is 3 1. 4 kms" 1 . This is c l ose to the values of 



v n t± < 15 - 30 kms 1 found by 



Hahn et al. 



(|2010| ). lLandi fc Cranmerl (120091 ) performed a 



similar analysis, but with the condition T- h ±_ > instead of T- h ± > T e , and found t> n t,_i_ < 30 
- 35 kms" 1 . Using the same conditions as theirs for T^± to reanalyze our data, we find 
Vnt, i < 28.3 kms" 1 . Thi s scales to 35.7 kms" 1 at 1.05 R Q , in reasonable agreement with 



Landi fc Cranmerl (120091 ) . Thus, the values for i>nt,± found on-disk appear consistent with 



the off-limb observations. 



3.3. Parallel Broadening 

The relation among f e ff,|l> and v at u is analogous to equation (Q. The thermal 
and non-thermal contributions in this direction can be estimated using the fact that 
both stoc hastic or ion cyclotron resonant heati ng are predicted to be weak in the parallel 



direction (ICranmer et al. 



1999 



Chandran 



20101 ). and so is set by collisions with protons 



and electrons. 



In the low corona electron-ion collisions are common. The Spitzer electron-ion 



- 19 - 



temperature equilibration time is ( ISpitzerl Il962l ; iBochslerl 120071 ) 



3Mm P 



eq ' e 327r 1 /% e Z i 2 e 4 lnA V M 



2k B Ti 2k B T e 



3/2 



(7) 



where Z\ is the ion charge, e is the elementary charge, m e is the electron mass, and In A 
is the Coulomb logarithm, which is ~ 20 for the corona. To estimate r eqie we can take 
n c = 2 x 10 8 cm~ 3 , and 2] = T e = 10 6 K. The electron-ion temperature equilibration time 
is then r eq e R3 30 s. Because of their more similar masses proton-ion collisions result in 
a more efficient energy exchange so that the proton-ion equilibration time T eqp is about 
^m c /m v shorter than r cq>e , or about .7 s here. This ca n be compared to the outflow 



timescale through the observed height. 



Tian et al. 



( I2010l ) has found the outflow speed in 



on-disk coronal hole measurements to be < 10 kms -1 . We estimate the height range of the 
observation to be > 0.01 R & . Thus, the outflow timescale is r out > 700 s, which is much 
longer than the equilibration time. This argument is supported by computational models 
for ion heating, which a re consistent with Ty = T p and isotropic T p a t such low heights in 



the corona of < 1.3 R Q ( jCranmer et al 



1999 



Landi Sz Cranmer 



Since parallel heating is predicted to be weak, we expect all the ions should be in 
equilibrium with the protons and all the ions should have the same 7] n = T p . Additionally, 
if the emission comes from the same volume then all the ions are subject to the same fluid 
motions so that v n t,\\ is the same for all the ions. To test whether these conditions are 
consistent with our data we found values for t> n t,|| and Ty that match the measured tyy 
for each measured ion species to within 2a or better. This test shows that v nt) \\ = 6.0 - 
14.9 kms -1 and Ty = 1.5 x 10 6 - 2.1 x 10 6 K are consistent with the data. Note that, as in 
the perpendicular case, for the analysis we implicitly assumed that the possible presence of 
multiple structures within each bin do not affect the analysis. This assumption is justified a 
posteriori by the consistency of our results. A possible reason that the higher T e structures 
do not affect the analysis is that their v nt 11 and Ty may not be too different from those of 
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the cooler structures. This could be caused by having different electron heating, but similar 
waves and proton and ion heating in the different structures. 

Having established that nearly uniform values for the ion temperature and nonthermal 
velocity are consistent with the data, we determined the most probable values by performing 
a least squares fit to the parallel analog of Equation [2] using the measured f e fr,||- In 
this fit the known ion masses M are the only independent variables. We find that 
Ty = (1.8 ± 0.2) x 10 6 K and v nt; \\ = 12.6 ± 2.3 kms" 1 . 

Figure [7] shows Ty versus T t after setting v ntt \\ to the fitted value of 12.6 kins" 1 . The 
dotted line indicates the fitted Ty and the solid line illustrates Ty = T t . The plot shows 
that the inferred Ty is greater than T t over a wide range of formation temperatures. 

One possible interpretation for T; n > T t is that the ions reflect a proton temperature 
that is roughly twice T e . Our Ty value does fall within the range of estimated proton 
temperatures in the corona, though there are few measurements of T p at low heights. 



Marsch et al. 



(120001 ) have measured T p « 2 x 10 5 K at ~ 1.02 R Q based on hydrogen line 
widths. We have estimated the height of this observation to be within about 0.03 R® 
above 1.02 R & . Thus, our data imply a large temperature gradient in this height range. 
Measurements with UVCS of the Lya line above 1.3 R^ hav e shown T p > T e , with a value 



of up to 3 x 10 6 K flEsser et al 



1999 



Antonucci et al 



20001 ). Our estimate for Tj ii falls 



between the proton temperature measurements at slightly smaller and somewhat larger 
heights. This may provide a useful constraint for solar wind models. 

Comparing Ty to T^±, we find that because of the large spread in the lower and 
upper bounds for T^± the ion temperatures in the parallel and perpendicular directions are 
similar. This may be due to either strong collisional cooling of the ions, most likely by 
protons, or weak perpendicular heating at these heights. Only for small q/M < 0.18 are the 
data consistent with T{± > 7] n. In contrast, UVCS measurements at 2.2 R Q have shown 
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T; 5 j_ > Tjii for O VI (q/M = 0.31). This implies that perpendicular heating of higher q/M 
ions becomes more efficient at larger heights. 

The parallel non-thermal velocity is th ought to be due to slow magnetosonic waves 



( jMcLements et al. 



1991 



Erdelyi et al. 



19981 ) . In addition to broadening the line widths, 



slow waves a re predicted to produc e asymmetric line shapes, which we did not find in this 



observation. 



Verwichte et al 



( 120101 ) has, however, shown that for low amplitude waves with 
v m,\\ 15 kms -1 the asymmetry would be only a few percent. Thus, for our inferred i> nt) n 
of 12.6 ± 2.3 kms -1 the asymmetry is predicted to be small, which may explain why we did 
not detect any. 

Magnetosonic waves are com pressional, and can t herefore be observed as time 



varying line intensity oscil 



ations (IBanerjee et al. 



20111 ) . In polar coronal hole plumes, 



DeForest fc Gurmanl (119981 ) found intensity variations corresponding to density amplitudes 
of 5% - 10% propagating n ear the sound speed. Similar perturbations in interplume 



regions have been found by 



Banerjee et al 



(l200lh with amplitudes also of a few percent. 



These results suggest that the observed oscillations are due to slow magnetosonic waves. 



Banerjee et al. 



( 1200 ll ) estimated the amplitude for the observed waves to be ~ 10 kms x . 



This amplitude is consistent with our inferred value for v nt \\. 



3.4. Outflow Velocity 

Combining measured Doppler shifts with the known magnetic field direction provides 
a new method for determining the outflow velocity. The velocity projection along the line 
of sight of the bulk flows parallel, v\\, and perpendicular, v±, to the magnetic field is given 
by ^los = v\\ cos a + v± sin a. This line of sight velocity produces a Doppler shift given 
by Aq — A = Ao(fLos/ c )) where Aq is the rest wavelength. Note that with this definition a 
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positive velocity corresponds to a blueshift. For the corona, the transverse motion of the 
magnetic field lines is expected to be small, v± ~ 0. Thus, A and vn can be determined 
from a linear least squares fit to A versus cos a: 



A = Ao — Ao — cos a. 

c 



(8) 



An example fit is shown in Figure |HJ Differential rotation could cause systematic errors in 
the analysis since a varies with latitude. However, based on the rotation rates measured by 



WohletaL 



(120101 ) we estimate that the rotation velocity differs by less than 0.2 kms -1 over 
the latitude range in our data. We therefore neglect differential rotation in the analysis. 
The rest wavelength on the detector Ao and its la uncertainty were found directly from 
the least squares fit, and these values were used in deriving vn and its uncertainty from 
the fitted slope. We fou nd that, on average, the inferred A n agreed with the wavelengths 



tabulated in CH1ANTI flDere et al 



1997 



Landi et al. 



2012ah to within 0.004 A. 



Figure [9] plots the inferred vn as a function of the DEM-averaged formation temperature 
logTi. For the analysis we performed the fit to equation (JHj) for each line separately to 
determine v» and then within each ion species took the weighted mean and error. The data 
point at logT t = 5.26 comes from the O IV line, which is not used in the line width analysis. 
Since these lines come from the transition region it could be that a two-Gaussian fit would 
make a better fit, however the estimated systemati c error from using a single Gaussian to 



fit these lines is ~ 0.5 kms 1 ( Peter fc Judg< 



19991 ) 



There is a systematic uncertainty in v» due to the variation of the detector wavelength 
scale as a function of the Solar-Y position. This is caused by the tilt of the EIS slits relati ve 



to the CCD. We corrected for this variation using the parameters given by 



Kamio et al. 



(120101 ) . However, due to the uncertainty in the fitting parameters, t here is a possible la 



systematic offset in the derived vn of about 1.3 kms 1 (jYoung 



2011M ). 



The results show that at coronal hole temperatures of logT e w 5.8 - 6.1 there is 
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an outflow velocity of w 5 kms^ 1 . The flow velocity is small at both lower and higher 
temperatures. The small velocity at low temperatures implies that the outflow is smaller 
in the transition region. There may be two reasons for this. First, in the transition 
region there could be low lying cool closed loops which have no net outflow. Second, if 
temperature is proportional to height, then the increasing v\\ with temperature may reflect 
acceleration within the transition region. The low velocity at high temperatures could be 
another indication that there are multiple structures within the large spatial bins and that 
the warm structures do not have strong outflows. 



£ ( 


1999 


) and 


Tian et al. 



( 1201 Of ). iPeter fc Judge! (119991 ) measured Doppler velocities at disk center and found 
outflow velocities of ~ 5 kins" 1 for logT e > 6.0, in reasonable agreement with the present 
measurements. However, they also found red-shifts at lower temperatures, in contrast to 
our results which show no flows at low temperatures. The reason for the di fference could b e 



that w e are observing open field lines while the disk center observations of 



Peter fc Judge 



(119991 ) are dominated by closed field lines. Material on the open field lines is free to escape, 



while the trapped material on t 
of a mass cycle as described by 



re closed field lines ma y return to the chromosphere as part 



Mcintosh et al. 



(120121 ). 



Tian et al 



( 120101 ) measured Doppler shifts from a polar coronal hole and found that 
the outflow velocity increased steadily with formation temperature, reaching ~ 25 kms" 1 
at temperatures of logT e = 6.2. In our analysis we find much smaller velocities that are 
constant over a b road temperature range. One possible reason for the discrepancy is that 



Tian et al. 



(120 101 ) focussed specifically on regions with a strong outflow. In the present 
analysis we performed no such selection. Since our method involves fitting over a range of a 
it necessarily finds an average outflow velocity for the coronal hole and could miss isolated 



regions of large outflow. The broad binning might also reduce the inferred velocity if such 
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vel ocities were confin ed to small regions. Another possible reason is that the observation 
by 



Tian et al. 



( 120101 ) looked near the boundary of the coronal hole, where there might be 
factors that increase the outflow, such as reconnection between the open and closed field 
lines. 

Our measured velocities are also much smaller than the velocities of pro pagating; 



Mcintosh! (J2012J). They 



coronal disturbances found in an equatorial coronal hole by 
used sequences of Atmospheric Imaging Assembly (AIA) images to identify propagating 
disturbances moving with apparent velocities of ~ 75 kins" 1 . Such intermittent events 
would appear as line asymmetries in our data, but as discussed above (Section I2.2p we did 
not find any significant line asymmetries. If these jets occur relatively infrequently or cover 
small spatial scales, then they would make up only a small percentage of the total emission. 
It is possible that we do not resolve the asymmetries due to insufficient signal to noise or 
the large spatial integration in our data. 



4. Conclusion 

We have described a new technique for combining magnetic field and spectroscopic 
data in order to derive plasma and wave properties in a coronal hole. We used a PFSS 
magnetic field model to determine the inclination angle between the line of sight and the 
magnetic field line in each spatial bin of our EIS spectrum. From these data we were able 
to infer the anisotropy of the ion temperatures, the non-thermal velocity, and the Doppler 
shift of plasma flowing along the field. The method provides both new results and also 
confirms other results that have been obtained using different methods. 

From the anisotropic line widths we inferred ion temperatures and non-thermal 
velocities induced by both transverse waves (Alfven) and longitudinal waves (slow 
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magnetosonic). The extracted t> n t,_i_ was consistent with estimates from off- limb observations 
for the Alfven wave amplitude and v nt » was consistent with observations of quasi-periodic 
intensity oscillations seen in coronal holes. The inferred Ti ± showed charge to mass ratio 
dependent heating, but with evidence suggesting greater collisional cooling than is found 
in off-limb measurements. The parallel ion temperature is expected to be equal to the 
proton temperature, for which there is limited data in the low corona. A comparison of 
our results for Ty with proton temperature measurements at slightly lower heights than 
this observation implies a large proton temperature gradient around 1.02 R & . Our Doppler 
shift results indicate an average outflow velocity of about 5 kms" 1 in the coronal hole over 
a broad temperature range. 

One limitation for the current data set was the large spatial binning that was required 
in order to obtain good statistics for fitting the line profiles. Our analysis shows evidence 
that some of the material in these bins comes from warmer than average structures within 
the coronal hole. An observation with higher spatial resolution and better statistics 
may resolve some of these systematic issues. Similar techniques could also be applied to 
structures o ther than coronal holes. For exa mple, nanoflare models predict anisotropic line 



broadening (IPatsourakos fc Klimchuk 



20061 ) which might be observed in active regions or 
coronal loops. Such an analysis would likely require a more sophisticated magnetic field 
model than a PFSS. 



We thank Enrico Landi for helpful suggestions and comments on the manuscript. 
MH and DWS were supported in part by the NASA Solar Heliospheric Physics program 
grant NNX09AB25G and the NSF Division of Atmospheric and Geospace Sciences SHINE 
program grant AGS- 1060 194. 




Fig. 1. — The box outlines the position of the EIS observation overlayed on a nearly contem- 
poraneous EIT SOHO image in the 195 A band, which consists mainly of Fe XII emission. 
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Fig. 2. — Fe XII 195.12 A line intensity measured by EIS. The grid shows the binning used 
in the analysis and the dots indicate those bins that were used in the analysis. These were 
chosen to avoid bright points and the quiet Sun regions in the bottom and rightmost portions 
of the observation. 
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Fig. 3. — PFSS model tracing the magnetic field lines that pass through the bins used in 
the analysis (Figure [2]). The image has been rotated 50° about the North-South axis in 
order to emphasize the changing angle of the field lines with respect to the line of sight. 
The dark curve indicates the surface at 1.05 R®. The grayscale on the solar disk shows the 
magnetogram data us ed for the model with whi te and black indicating positive and negative 



polarity, respectively ( ISchrijver &: Derosa 



20031 ) 
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Fig. 4. — Effective velocity v e $ versus the angle a between the line of sight and the magnetic 
field for the Fe XI lines. The solid line shows the fit to equation (J3]). 
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Fig. 5. — Representa t ive D EM for the coronal hole. The solid line shows the results from 



the 



Landi fc Landini ( 



19971 ) inversion method and the crosses show the results from the 



Hannah fc Kontarl (120121 ) method. The filled circles show the points on the DEM determined 



from the individu al measured line i ntens ities and their scatter gives an estimate of the 



uncertainty in the lLandi fc Landinil ( 1l997l ) inferred DEM. On the temperature scale, these 
points correspond to log T t for each line (see text). The size of the various error bars indicates 
the la uncertainties. 
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Fig. 6. — Upper and lower bounds for log 7]^ as a function of charge to mass ratio q/M. 
The top axis labels the corresponding ion species. 
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Fig. 7. — Parallel ion temperature log 7] ii as a function of the DEM-averaged formation 
temperature logT t . The filled circles indicate after setting v nt t \\ to the fitted value of 
f nt> || = 12.6 ± 2.3 kms" 1 . The dotted line shows the best fit temperature of = (1.8 ± 
0.2) x 10 6 K. The solid line illustrates where = T t . 
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Fig. 8. — Linear fit of A versus cos(o;) for the Fe IX 189.94 A line. For this particular line 
the fit finds A = 189.9362 ± 0.0002 A and v\\ = 5.9 ± 0.6 kms" 1 . 
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Fig. 9. — Inferred flow velocity v» along the magnetic field as a function of the DEM-averaged 
formation temperature logT t . The dotted horizontal line shows v\\ = 0. There is a possible 
systematic offset in v\\ of ±1.3 kms" 1 due to the uncertainty in correcting for the tilt of the 
EIS slits relative to the CCD (see text for more detail). 
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Table 1:: Line List. 





Ion 


A (A) 1 
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He II 


f 256.317 


ls 2 S 1/2 




3» 2 P 3/2 

o/Z 












ZdU.OlO 
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9 2 9 2 p 
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IVlg V 


97fi X.7Q 

z i o.o ( y 
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2s 2o 5 1 P 
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9 c2 o„3 2 n 
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O J 9„3 2 n 
ZS Zp JJ 3 /2 




2s2» 42 P,/o 
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97R 1 KA 
Z 1 0. 1041 


9„2 9rl 2 3 p 
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3s 2 3p 34 S 3/2 




3s 2 3p 2 ( 3 P) 3d 4 P 1/2 




Fe XII 


193.509 


3s 2 3p 34 S 3/2 




os op .r J oa / 3 / 2 


* 


Fe XII 


j- 195.119 


3s 2 3p 3 4 ,S 3/2 




3s 2 3p 2 ( 3 P) 3d 4 P 5/2 






1 195.179 


3s 2 3p 3 2 P 3/2 


— 


3s 2 3p 2 ( l D)3d 2 D 3/2 


* 


Fe xill 


202.044 


3s 2 3p 2 3 P 




3s 2 3p3rf 3 Pi 






203.772 


3s 3p 3 3 Pi 




3s,3p 2 3d 3 P 2 




Fe xill < 


203.796 


3s 2 3p 2 3 P 2 




3s 2 3p3d 3 D 2 






203.827 


3s 2 3p 2 3 P 2 




3s 2 3p3d 3 D 3 






203.835 


3s 3p 3 3 P 2 




3s,3p 2 3d 3 F 2 




Fe xiv 


270.521 


3s 2 3p 2 P 3/2 




3s3p 22 P 1/2 



-37- 



Ion 



A (A) 1 



Transition 1 



Fe XV 



284.163 



3s 2 l S n 



3s 3p 1 P 1 



1997 



LandietaL 



2012ah . 



1 Wavelengths and transitions taken from CHIANTI (IDere et al. 
°Brackets indicate blends from the same ion. Asterisks mark the lines used in the analysis 
of the line widths. These and the remaining lines were all used for the DEM analysis. 
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Table 2: Line Widths and Effective Formation Temperatures 



Ion 


v eSi \\ (kms -1 ) 


v eS>± (kms -1 ) 


logpt(K)] 


Mg vi 


36. f ± 1.7 


33.6 ±1.0 


5.68 


Mg vii 


43.2 ±5.4 


31.4 ± 1.8 


5.81 


Si vi 


33.8 ±8.0 


40.0 ±3.7 


5.67 


Si vii 


35. f ±0.8 


36.2 ±0.3 


5.82 


Si x 


30.7 ±2.5 


36.6 ± 1.2 


6.10 


S x 


42.6 ±6.4 


35.0 ±2.8 


6.12 


Fe vm 


26.5 ±0.7 


32.2 ±0.3 


5.75 


Fe IX 


25.5 ± f.l 


34.1 ±0.4 


5.91 


Fe x 


27.2 ±1.0 


32.4 ±0.4 


6.01 


Fe XI 


24.2 ±0.9 


29.1 ±0.4 


6.08 


Fe XII 


27.1 ±0.7 


30.7 ±0.3 


6.14 


Fe xill 


23.2 ±2.3 


27.5 ±1.0 


6.20 



-39 - 
REFERENCES 

Antonucci, E., Dodero, M. A., & Giordano, S. 2000, Sol. Phys., 197, 115 

Banerjee, D., Gupta, G., & Teriaca, L. 2011, Space Sci. Rev., 158, 267 

Banerjee, D., O'Shea, E., Doyle, J. G., & Goossens, M. 2001, A&A, 377, 691 

Banerjee, D., Teriaca, L., Doyle, J. G., & Wilhelm, K. 1998, A&A, 339, 208 

Belcher, J. W. & Davis, L. 1971, J. Geophys. Res., 76, 3534 

Bemporad, A. & Abbo, L. 2012, ApJ, 751, 110 

Bochsler, P. 2007, Astron. Astrophys. Rev., 14, 1 

Chandran, B. D. G. 2010, ApJ, 720, 548 

Cranmer, S. R. 2002, Space Sci. Rev., 101, 229 

Cranmer, S. R. et al. 1999, ApJ, 511, 481 

Culhane, J. L. et al. 2007, Sol. Phys., 243, 19 

De Pontieu, B. et al. 2007, Science, 318, 1574 

DeForest, C. E. & Gurman, J. B. 1998, ApJ, 501, L217 

Delaboudiniere, J. P. et al. 1995, Sol. Phys., 162, 291 

Dere, K. P., Landi, E., Mason, H. E., Fossi, B. C. M., & Young, P. R. 1997, A&AS, 125, 149 
Doyle, J. G., Banerjee, D., & Perez, M. E. 1998, Sol. Phys., 181, 91 
Erdelyi, R., Doyle, J. G., Perez, M. E., & Wilhelm, K. 1998, A&A, 337, 287 
Esser, R. et al. 1999, ApJ, 510, 63 



-40- 

Feldman, U., Doschek, G. A., & Patterson, N. P. 1976, ApJ, 209, 270 
Hahn, M., Bryans, P., Landi, E., Miralles, M. P., & Savin, D. W. 2010, ApJ, 725, 774 
Hahn, M., Landi, E., & Savin, D. W. 2011, ApJ, 736, 101 
-. 2012, ApJ, 753, 36 

Hannah, I. G. & Kontar, E. P. 2012, A&A, 539, 146 

Hara, H., Watanabe, T., Harra, L. K., Culhane, J. L., & Young, P. R. 2011, ApJ, 741, 107 
Hoeksema, J. T., Wilcox, J. M., & Scherrer, P. H. 1983, JGR, 88, 9910 
Hollweg, J. V. 2008, J. Astrophys. Astr., 29, 217 

Kamio, S., Hara, H., Watanabe, T., Fredvik, T., & Hansteen, V. H. 2010, Sol. Phys., 266, 
209 

Kjeldseth Moe, O. k Nicolas, K. R. 1977, ApJ, 211, 579 

Kohl, J. L. et al. 1995, Sol. Phys., 162, 313 

-. 1998, ApJ, 501, L127 

Landi, E. 2012, private communication 

Landi, E. k Cranmer, S. R. 2009, ApJ, 691, 794 

Landi, E., Del Zanna, G., Young, P. R., Dere, K. P., k Mason, H. E. 2012a, ApJ, 744, 99 
Landi, E. k Landini, M. 1997, A&A, 327, 1230 
Landi, E., Reale, F., k Testa, P. 2012b, A&A, 538, 111 
Marsch, E., Ao, X.-Z., & Tu, C.-Y. 2004, JGR, 109, 04102 



-41 - 

Marsch, E., Tu, C.-Y., & Wilhelm, K. 2000, A&A, 359, 381 
Mcintosh, S. W. 2012, Space Sci. Rev., in press 
Mcintosh, S. W. & De Pontieu, B. 2009, ApJ, 707, 524 

Mcintosh, S. W., Tian, H., Sechler, M., & De Pontieu, B. 2012, ApJ, 749, 60 

McLements, K. G., Harrison, R. A., & Alexander, D. 1991, Sol. Phys., 131, 41 

Moran, T. G. 2001, A&A, 374, L9 

Neugebauer, M. et al. 1998, JGR, 103, 14587 

Patsourakos, S. & Klimchuk, J. A. 2006, ApJ, 647, 1452 

Peter, H. 2001, A&A, 374, 1108 

Peter, H. & Judge, P. G. 1999, ApJ, 522, 1148 

Phillips, K. J. H., Feldman, U., & Landi, E. 2008, Ultraviolet and X-ray Spectroscopy of 
the Solar Atmosphere (Cambridge University Press) 

Riley, P., Linker, J. A., Mikic, Z., Lionello, R., Ledvina, S. A., & Luhmann, J. G. 2006, 
ApJ, 653, 1510 

Roussel-Dupre, R., Francis, M. H., & Billings, D. E. 1979, MNRAS, 187, 9 

Schatten, K. H., Wilcox, J. M., & Ness, N. F. 1969, Sol. Phys., 6, 442 

Schrijver, C. J. & Derosa, M. L. 2003, Sol. Phys., 212, 165 

Spitzer, L. 1962, Physics of Fully Ionized Gases (John Wiley and Sons, Inc.) 

Tian, H., Mcintosh, S. W., De Pontieu, B., Martmez-Sykora, J., Sechler, M., & Wang, X. 
2011, ApJ, 738, 18 



Tian, H., Tu, C, Marsch, E., He, J., & Kamio, S. 2010, ApJ, 709, L88 
Tomczyk, S. et al. 2007, Science, 317, 1192 

Tu, C.-Y., Marsch, E., Wilhelm, K., & Curdt, W. 1998, ApJ, 503, 475 

Verwichte, E., Marsh, M., Foullon, C, Van Doorsselaere, T., De Moortel, I., Hood, A. W., 
& Nakariakov, V. M. 2010, ApJ, 724, L194 

Voitenko, Y. & Goossens, M. 2004, ApJ, 605, L149 

Wang, Y.-M. & Sheeley Jr., N. R. 1992, ApJ, 392, 310 

Wilhelm, K., Marsch, E., Dwivedi, B. N., & Feldman, U. 2007, Space Sci. Rev., 133, 103 

Wohl, H., Brajsa, R., Hanslmeier, A., k Gissot, S. F. 2010, A&A, 520, 29 

Young, P. 2011a, EIS Software Note No. 7: 

http:/ /hesperia. gsfc.nasa.gov/ssw/hinode/eis/doc/eis_notes/ 
eis_notes/07_LINE_WIDTH/eis_swnote_07.pdf 

-. 2011b, EIS Software Note No. 4: 

http : / /hesperia. gsfc . nasa. gov / ssw/hinode /eis/doc/ eis_notes / 
eis_notes/04_SLIT_TILT/eis_swnote_04.pdf 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



